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INTFODUCTiaJ 


This report concems the rotordynairdos of the HPFTP (High Pressure Fuel TU3±opuitp) 
of the SSME (Space Shuttle Main Engine ) , and is a sequel to several earlier studies of 
this unit by the author [l], [2], [3]. The results of both linear (stability and 
synchronous response) and transient nonlinear analyses are reported. The same anal^^is 
procedures are enroloyed here as previously; however, the data defining the rotordynanicJ 
model differs in the following two regards: 

(a) The most-recent MSFC case structural dynamics model is enployed. 

(b) New dynamic coefficients were developed in the course of this study for tlie 
HPFTP interstage seals [4], [5], and introduced into the xotordynamics model. 
A copy of reference [5] is enclosed as part of this final report. 

The objective of the analysis reported herein was an examination of possible 
problems in moving frcrn KPL to FPL conditions, and a review of prior analysis results 
with tlie 'dpdated rotordynamics model as described above. The following qi^tions are 
examineci by the analysis corpleted: 

(a) What would be the influence on HPFTP rotordynamics of a change in inter- 
stage seals from the presently orployed "smooth-stepped" design to a 
"sanoothrstraight" configuration? 

(b) How sensitive are stability and synchronous results to changes in bearing 
stiffnesses and danping? 

(c) What would be the influence on rotorc^maraic stability of a change frcm the 
present stiff synmetric bearing-carrier design to an asyrrmetric bearing- 
carrier configuration? 
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TOE HPFTP ROTORDYNAMIC ,®DEIi 


The rotordynainic nodal used in this stu(^ basically resonbles those errployed by 
the auti’sor [l], [2], [3] and other rotordynamics investigators of the SSME turbopurrps. 
The data and paraneters vAiidi aire required to define a rotordynamics model can be 
separated into those vdiich are relatively well known and assumed fixed, and those 
which are known only within limits, and are to be varied. 

Specified Dat a 

In the present analysis, the follodng parameters are assumed to be known and 
fixed; 

(a) rotor and case structural dynamic models, 

(b) speed-dependent dynamic coefficients for seals, 

(c) stiffness and damping coefficients for the balance piston, 

(d) bearing carrier stiffnesses, 

(e) bjaririg dead-band clearances, and 

(f) ii±»alance magnitude and location. 

In a recent study of the HPOTP (High Pressure Oxygen Turbopump) , the bearing 
dead-hand clearances coiiLd 'vary widely due to tolerance stackups, This is not the 

case with the HPFTP which operates with much snaller clearances on the order of 

—3 

0.25 X 10 inches. A 6 gm-in imbalance magnitude, located between tha turbine 
wheels, was used to represent the effect of irregular coating loss on the turbine 
blades. Appendix A provides the numbers used to define the balance of the data. 

Hie sources for the data are as follows: The rotor structural dynamic model was 

developed by B. Rowan, Rocketd^Tie; while the case structural dynamiic model was developed 
by MSEC. As noted earlier, the s^l coefficients were developed by the'.^giithor as 
documented in [4], [5]. Hie ranaining data coincide with those currently employed by 
MSFC as provided by S. Winder in an annotated cxxrputer data listing. 

Varied Data 

Hie ncndnal radial stiffness model of the HPFTP bearings as a function of speed 
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was provided by S. Winder, and is bas«3 on A.B. Jones analysis results as nodified 
due to e:xp^inental data. These bearing stiffnesses seem to be reasonable in that 
they yield the correct first-critical-speed location. However, the PIPFTP bearings 
are the sane as the HPCXCP (High Pressure Oxygen Turbopunp) prebumer bearings, and 
estimates of the SPOTP bearings I'iave varied markedly. Hence, in this stucty, the 
stiffness is varied by +25% of nctninal. 

Bearing damping is modeled by a linear, constant, danping coefficient with the 
values Cg=0,20 lb sec/ro. 

A "custonary" destabilizing mechanism acting on a rotor's turbines is the "Alford 
aerodynamic cross-co\:pling" force [6] which is modeled by 
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viiere T is the turbine torque, is the average pitch diameter of the turbine blades, 
H is tf’je average height, and is the "change of thermodyanmic efficiency per unit of 
rotor displacement, expressed as a frac±ion of blade height." Ihe physical rationale 
for tlieye forces is based on an increase in blade efficiency wi^di decreasing tip 
clearance. The parameter $ is varied here from uniiy to three. 

No attenpt is made to model iirpeller-diffuser forces similar to those celcsulated 
by J. COlding- Jorgensen [7] for a plain volute configuration, and mere recently for 
the vaned-dif fuser of the HPOTP. No such (calculations have been made for the HPFTP 
irtpeTilers. Since these forces are assxmed to be proportional to density, they will 
clearly be smaller for the HPFTP than the HPOTP. 

Based on earlier results mth the HPOTP [8], the influence of drive and load 
torques is rot included in the present study. 
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KENEMl STmiLTFi m> SYNCHRONOUS BESPONSE RESULTS 


Introduction 

This section provideij results based on linear synchronous response and stability 
calculations enploying the procedures outlined in reference [2]„ TtfS results 
presented should be interpreted in li^t of the kncwn nonlinear effects of bearing 
clearances [9]. Specifically, bearing clearances have a softening effect in lowers 
ing the apparent critical speed location, and cause a reduction in the peak magnitude 
of rotor response at the nonlinear* critical speed, as opposed to the zero- 
clearance linear critical speed. Ifciwever, at ninning speeds below the nonli ne ar 
critical speed, airplitudes are increased with increasing bearing clearances. Note 
further in examining the linear results of this section that the axial balance- 
piston coiipling between the rotor and case is not included. However, this is included 
in the transient nodel, viiose results are included in the next 

Figure 1 i3.1ustrates the reference coordinate system. The Z ^ds of this 
figure is seen to be aligned vdtii the nominal axis of rotation for the HPFTP. As 
explained in [2], stability and synchronous re^nse calculations are based on 
coupled rotor-case .mades, which are calculated frcKi separate case data (excluding 
the HPPTP rotor) , fr^-free rotor modal data, and limar bearing stiffnesses. 

Sepeirate modes are calculated for both the X-Z and Y-Z planss* 

Orthotropic Bearing Carrier Design 

Previous rotort^Tiamics analysis of the HPPIP suggest that stability could be 

markedly iitproved by using an ortliotropic bearing carrier with orthogonal 

5 c 

stiffress(!ss of 1.75 X 10 and 5.0 X 10 Ibs/in respectively. This configuration 
requires bearing damping, and the original grooved interstage seal design. 

Listed below are the orthotropic support and daitping cases examined, together with 
their onset speed of instability results. 
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*Peak-artplitude running speed. 





Figure 1. Loce^ coordinate system for HPFTP analysis 
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(a) K =1.75 X 10^ Win, K =5.0 X 10^ Ib/iji 
C =5 Ib sec/in: w =26,110 xpcti 


V20 lb *=/“= 


w =49,710 xpn 
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(b) K =5.0 X 10^ Win, K =1.75 X 10^ Ib/in 
C =5 lb sec/in: o) =26,140 rpm 


C «20 lb sec/in: 
s 


( 1 ) =49,170 rniv 
s 


Ihe c3anpdj]ig support constant C corresponds to a linear "dash-pot" in parallel 
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vdth the sr^jport stiffness. 

Ihese stability results are not as attractive as previous predictions for 


ortli'irotrcpic bearing supports. Specifically, in pri.or analyses siroport orthotropy 
caused the rotor to be stable relatively independent of ivi^port dartpii^. With the 
present relatively-soft structural dynamic model, rotordynaraic stabiliiy is primarily 
acliisved by .increased daitping, and not by orthotropy. Note that stability results 
are relatively indifferent to the alignment of the bearing orthotropy axes with 
either the X-Z or Y-Z planes. 

Figures 2 and 3 illustrate synchronous response results for C =20 lb sec/in. 
Bearing reactions are seen to be quite low for either alignment of the bearing 
orthotropy a:<es. In these figures, and throughout this report, bearing nunbers 1 
through 4 refer, respectively, to the punp outboard, punp inboard, turbine inboard, 


and turbine outboard bearings. 

Conventional, Symmetric, Bearing-Carrier Design 

Table 1 summarises stability results for stepped and straight interstage seal 

designs for a 25% variation in noninal bearing stiffness with the two support dairp- 

ing coefficients, C =5,20 lb sec/in. The entry C(FPL) denotes the percent of 

s 

critical damping at FPL, with w continuing to denote tte onset speecs of instability,. 

s 

The results pcresented are insensitive to changes in support damping. Except 
for the soft 0.75 Kj^ / stepped -seal combination, the onset speed of instability 
results are largely insensitive to changes in bearing stiffness or interstage seal 
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(a) bearing reactions (lbs) 

(b) punp accel. levels (g*s) 

(c) turbine acc^l. levels (g's) 
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stepped Seal 


0.75 

Cjs5 

Wg«74,544 

C (FPL) *1.15% 



uj^=*75#OOQ 

15 

5 (FH.)=1.16% 

1.0 


w =54,733 

5 

C (PPL) =0.88% 


C =20 
“ 

ti) =54,900 

5 

C (FPL) ^.89% 

1.25 K 

b 

C «5 
s 

o> =53,663 
s 

C (FPL) =2. 9% 


C-20 

s 

Og=54,S97 

C (PPL) =2. 5% 


Straji^t Seal 


0.75 

"b 


03^=51,634 

5 (FPL) =1.5% 



0^=20 

0) =51,890 
s 

? (FPL) =1.5% 

1.0 



03 =52,605 
s 

C (PPL) =1.5% 



C =20 
s 

03 =52,833 
s 

C (FPL) =1.6% 

1.25 

*'b 

C=5 

03=53,782 

s 

e (FPL) =1.6% 



C=20 

03 =53,930 

C (FPL) =1.6% 


Table 1. Onset speed of instability 
results for the syntnetric bearing 
carrier vdth stepped and straight seals 
for 25% variation in bearing stiffeesses. 
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<3esigns, Uie interstage seal design becomes slightly irore stable with reduced 
bearing stiffness , while the converse situation holu» for the straight seal design. 

Figure 4 illustrates the reduction in onset sp©^ of instability for the stepped 
and interstage seal designs due to an increase in the Alford cross-coupling t^rm. 
These results are for tie nominal bearing stiffness coefficients C =5 lb sec/in. 

S 

Note that t>e straight-seal design is better able to withstand increases in 3 
than thn stepped seal, although either design has a ccxnfortable speed margin for 
stability. 

Figures 5 tteough 7 illustrate synchronous imbalance results for the stepped 
interstage seal design for 0.75 Kj^, 1.0 ajnd 1.25 bearing stiffness defini- 
tions. These results are for 6 girrin irrbalance between the turbine vheels, and 
oorrespcmd to support dartping of '0^-5 lb sec/in. The peak bearing loads are pre- 
dicted for a critical speed in the vicinity of 32,000 rpm, and are moderately 
reduced fox' the stiff bearing (1.25 Kj^) design. The purrp ir±>oard bearing is 
predicted to ex|»rienc8 the largest load on the order of 2000 lbs. However, the 
associated turbirie accel. levels are predicted to be betv^een 20 and 60 g's, vhich 
is markedly higher than measiarsd results for the HPFTP. 

Figures 8 through 10 provide the sarna results for the straight-seal configmra- 
tion as figures 5 through 7 provide for the stepped-seal. The bearing reactions 
loads are markedly reduced by the change in seals; hoover, the a<xel. levels are 
not proportionately reduced, and continue to be predicted at much higher levels 
than those which have been measured. The straight interstage seal configuration 
is seen to be quite insensitive to bearing stiffness variations. 

In examining the results of figures 5 thrcugh 10, recall that the bearing 
clearance effect is not included. The bearing dead-band will markedly reduce the 
peak bearing loads. The following conclusions may be re a c h ed fran an examination of 
these figures: 
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ONSET SPEED OF INSTABILITY (RPM) 


50.000 

45.000 

40.000 
FPL 

35.000 

30.000 
Figure 4 
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(a) bearing reactions (lbs) 

(b) puitp accel. levels (g’s) 

(c) turbine accel. levels (g's) 
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(b) ponp accel. levels (g's) 

(c) txarbine accel. levels (g's) 
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(a) bearing reactions (lbs) 

(b) piatp accel. levels (g’s) 

(c) turbine accel. levels (g®s) 



(a) l!he straight interstage seal configuration is preferrable to the stepped, 
because it yields iterkedly smaller bearing loads. 

(b) With the current case structural dynamics model, peak bearing loads are 
predicted at speeds around 32,000 rpm, vMdi is belcw I-<PL. Hence, the 
linear ^>n<^hronous response results actually predict smaller loads at 
FPL than FPL and lower speeds. 



TRANSIENT lOJLINEAR ANALYSIS RESULTS 


As noted in the introduction section, the principal concern in this study is 

potential problems to be encountered in moving fran RPL to FPL operating conditions. 

The analysis results reported in the preceding section address this concern via linear 

analysis procedures, vhile the results presented here consider the same question using 

transient nonlinear analysis techniques. The essential physical noiilinearity in the 

-3 

turbopunp is the bearing "dead-band" clearance on the order of 0.25 X 10 inches. 

Simulations were conducted for the nominal 1.0 bearing stiffnesses, zero 
bearing daitping, and straight and stepped-seal configurations. The results of a 
uniform deceleration from steac^-state initial conditions at FPL to RPL for the 
stepped and straight configurations are illustrated in figures 11 and 12, respectively. 
The stepped-seal configuration shows a sharp peak in the vicinity of 35,800 rpm vdth 
high associated bearing loads on the order of 2000 to 24CQ lbs; however, the 
associated punp and turbine accel. levels are not excessive, approaching 12 g's on 
the purtp 90 accel. output. The results presented in figure 12 for the straight-ssal- 
ooefficient simulation are much none encouraging. The critical in this case occurs 
at approximately 37,200 rpm, with ps'ak bearing loads on the order of 600 rpm. 

Peak bearing loads for the stepped seal configuration were calculated for a 
steady running speed of 35,808 rpm with the following results: 

R^ = 2000 lbs = 2300 lbs 

= 2200 lbs R^ = 2400 lbs 

The associated putip and turbine accel. levels for this steady-state operating condi- 
tions are: 

G(Punp rad.) = 3 g’s G(Turbine 90°) = 1.0 g's 

G{Punp 90°) -10 g's G(Turbine 180°)= 2.0 g's 

G(Punp 180°) =10 g's 
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Peak bearing loads for the straight-seal configuration were calculated at PPL 


conditions with the following results; 

= 400 lbs = 500 lbs 

R 2 = 400 lbs R^ = 500 lbs 

The associated punp and turbine accel. levels are: 

G(Punp rad.) = 2.0 g's G(Turbine 90°) = 1.2 g's 

G(Purtp 90°) = 7.0 g's G(Turbine 180°)= 1.5 g's 

G(Purrp 180°) = 6.0 g's 

Ihe unmistakeable conclusion to be reached fron these results is that predicted 
bearing loads are markedly reduced for the straight seal as opposed to the stepped- 
seal coefficients. The basic question \\tdch remains to be answered? however, is 
"Hew good are the calculated seal coefficients?" To the extent that they are an 
accurate measure of the seal's relative dynamic properties, the straight seal is 
obvioiisly to be preferred. 
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Figure 11 : Simolation of deceleration from FPL co RPL for 

the stepped -seal configuration. 
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Figure 12 : Sijmilation of deceleration fron FPL to RPL for 

the straight-seal configuration. 







Fi^e 12 ; Simulation of deceleration from FPL to FPL for 
the straight-seal configuration. 
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Figure 12 : Simulation of deceleration fran FPL to RPL for 

the straight-seed configuration. 
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APPENDIX A 


INPUT DATA FOR THE HPFTP PDODRDYNAMICS MODEL 

The fixed data used to define the rotordynamics model are provided in this 
appendix. Inch-lb-*sec units are used throughout. 

Rotor Eigenvalues 

TI,ie rotor eigenvalues and eigenvectors used here are based on a structural- 
dynamic model by B. Rowan. Hie free-free eigenvalues used are listed below: 



= 0 

^5 

= 2048.6 

Hz 

^2 

= 0 

CTl 

= 2622.7 

Hz 

^3 

= 632.72 Hz 

X7 

= 3155.4 

Hz 

^4 

=1397.2 Hz 

^8 

= 378**. 7 

Hz 


Case Eigenvalues and Dairping Factors 

The case eigenvalues and eigenvectors are based on the most recent MSFC 
structural dynamic model. The eigenvalues used in this study are: 


'cl 

= 271.04 

Hz 

^ c = 

c6 

561.79 

Hz 

'c2 

= 370.11 

Hz 

0 

II 

564.99 

Hz 

'c3 

= 440.28 

Hz 

0 

CO 

II 

609.84 

Hz 

'c4 

= 500,54 

Hz 

'c9 “ 

706.10 

Hz 

'c5 

= 512.59 

Hz 

'clo"" 

730.64 

Hz 


One-half percent of cri.tical danroing was used for all modes. These are the same 
modes used in the MSFC transient model. 
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Seal Dynamic Coefficients 


The only seals of inportance to the rotordynamic response of the HPPTP 
are the interstage seals. The coefficients for these seals are provided in 
Table 2(a) and Table 2(b) . 


K 

k 

C 

c 

M 


FPL 

5.003 X 10^ 
0.8187 X 10^ 
1D2.3 
13.34 

6.242 X 10“^ 


RPL 

4.3981 X 10^ 
0.7986 X 10^ 
94.2 
13.10 

6.033 X lO"^ 


MPL 

2.4242 X 10^ 
0.3275 X 10^ 
64.8 
7.562 

5.144 X 10“^ 


Table 2(a) . Dynamic seal coefficients for straight smooth seals. 


FPL 

K 1.247 X 10^ 

k 1.377 X lo"^ 

C 20.49 

c 1.317 

M 5.374 X la"^ 


RPL 

1,098 X 10^ 
1.191 X lo"^ 
18.91 
1.194 

5.197 X lO”"* 


MPL 

0.5044 X 10^ 
0.5458 X lO"^ 
12.99 
0.688 

4.939 X 10"^ 


Table 2(b) . Dynamic seal coefficients for smootli stepped seals. 

Balance-Piston Stiffness and Dynamic Coefficients 

The balance piston stiffness and dairping coefficients are: 

KZ = 2.6 X 10^ , CZ = 450.1 

Bearing & Bearing Carrier Stiffness 

The same bearing stiffness is used for all bearings, and is interpolated frcm 
the following data: 
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Speed 


Zero 


MPL 


RPL 


FPL 


''b 


1.38X10® 0.705X10® 0.654X10® 


The bearing support stiffness used is: 

Kg = 2.65 X 10® Ib/in. 


0.648X10 


DYNAMIC ANALYSIS. OF TURBULENT ANNULAR SEALS 
BASED ON HIRS' LUBRICATION EQUATION* 

D. W. Childs 

Mechanical Engineering Department 
Speed Scientific School 
The University of Louisville 
Louisville, KY 40208 

Expressions are derived which define dynamic coefficients for 
high-pressure annular seals typical of neck-ring and interstage 
seals employed in multi-stage centrifugal pumps. Completely 
developed turbulent flow is assumed in both the circumferential 
and axial directions, and is modeled in this analysis by Kirs' 
turbulent lubrication equations. Linear zeroth and first-order 
"short-bearing" perturbation solutions are developed by an expan- 
sion in the eccentricity ratio. The influence of inlet swirl is 
accounted for in the development of the circumferential flow field. 
Comparisons are made between the stiffness, damping, and inertia 
coefficients derived herein based on Kirs' model and previously 
published results based on other models. Finally, numerical 
results are presented for interstage seals in the Space Shuttle 
Main Engine High Pressure Fuel Turbopump. 


*This work was supported in part by NASA Grant NAS8-31233 from the 
George C. Marshall Space Flight Center and NASA Grant NSG 3200 
Lewis Research Center. 



Nomenclature 

a 

b 

B 

C/C 

C 

E 


h 



L 

m^ = 0.066, Hq = -0.25 
P 

Pi 

Po'Pi 

Ap = (l+C+2a)pvV2 
R 

= p (Rw) h/u 
R^ = pVh/y 
R^q = p(Rw)C/y 

U = RO) 

^ z '^0 

0-g = Ug/RC 


Dimensionless coefficient defined in 
Eq. (11) 

Dimensionless coefficient defined in 
Eq. (7) 

Dimensionless coefficient defined by 
Eq. (20) 

Dimensionless seal damping coefficients 
defined by Eq. (28) 

Nominal seal radial clearance (L) 

Dimensionless coefficient defined by 
Eq. (29) 

Seal radial clearance (L) 

First order perturbation in h (L) 

Drmensxonless seal stxffness coef f xcxenns 
defined by Eq. (28) 

Seal length (L) 

Coefficients for Hirs ' turbulence equations 

Fluid pressure (F/L^) 

Supply pressure at entrance (F/L^) 

Zeroth and first-order pressure 
perturbation (F/L^) 

Nominal pressure drop across seal (F/L^) 
Seal radius (L) 

Circumferential local Reynolds number 

Axial local Reynolds number 

Nominal circumferential Reynolds number 

Nominal axial Reynolds number 

(L/T) 

Axial and tangential fluid velocity 
components (L/T) 

Dimensionless velocity components 



“ao'“eo 


V 


V 

z /R6 
Z = z/L 

x,y 

e 

e 

C 

P 

X 

a = XL/C 
T = t/T 


zeroth order perturbetions in 
First order perturbations in U 

z d 

Dimensionless inlet swirl defined as 
the solution to Eg. (10) 

Initial (z=0) swirl 

Nominal axial velocity (L/T) 

Seal coordinates illustrated in 
figure 2 (L) 

Dimensionless axial seal position 

Radial seal displacements (L) 

Dimensionless coefficient defined by 
Eg. (11) 

Eccentricity ratio introduced in Eg. (1) 
Inlet pressure-loss coefficient 
Fluid density (ML“^) 

Friction loss coefficient defined in 
Eg. (11) 

Dimensionless time 

Shaft angular velocity (T^M 


tu 



Introduction 


In a series of publications, Black and Jenssen [1,2, 3,4] have 
explained the influence of seal forces on the rotordynaraic behavior 
of pumps. Figure 1 illustrates the two seal types v;hich have the 
potential for developing significant rotor forces. The neck or 
wear-ring seals are provided to reduce the back leakage flow along 
the front surface of the impeller face, while the interstage seal 
reduces the leakage from an impeller inlet back along the shaft to 
the back side of the preceding impeller. Black and Jenssen' s work 
dealt primarily with boiler-feed pumps handling water. However, 
seals may have a significant influence on rotordynamic behavior of 
any multi-stags pump, e.g. , reference [5] discusses the crucial 
dependence of the dynamic response of a hydrogen turbopump on its 
interstage seal design. 

Black and Jenssen are responsible for the initial analytical 
development of dynamic stiffness and damping coefficients for high- 
pressure annular seals. Their coefficients are valid for small 
motion of the seal journal about a centered position relative to 
its bearing. A bulk-flow analysis is employed, with the circum- 
ferential bulk-flow velocity assumed to be fully developed shear 
flow at The axial-flow momentum equation agrees with Yamada's 

[6] friction-loss results for rotating concentric cylinders, which 
defines the axial friction factor as a function of the axial and 
radial Reynolds numbe7.‘s. In [4], Black and Jenssen define the 
friction factor as a function of the local axial and radial Reynolds 
numbers . 

While Black and Jenssen 's results apply only for s;Tiall seal 
motion about a centered position, Allaire et al. [7] use Black's 



model to numerically calculate dynamic coefficients at large 
eccentricities. Further/ while Black and Jenssen define seal 
coefficients in a half-speed rotating reference system/ and employ 
a coordinate transformation to achieve stationary-reference results » 
Allaire et al. perform all calculations in a stationciry reference 
frame. 

In an as yet unpublished result [8] / Black has combined his 
prior seal-analysis governing equations with equations previously 
derived [93 for the analysis of, "Journal-bearings with high axial- 
flow in the turbulent regime/" to examine the development of cir- 
cumferential flow in a centered seal as a function of axial seal 
position. His results demonstrate that the bulk circumferential 
velocity of a fluid element asymptotically approaches -y- as it 
proceeds axially along the seal. The rate at which it approaches 
this limiting velocity depends on the amplitude of the axial and 
radial Reynolds numbers. Predictions of the stiffness cross- 
coupling term are generally reduced if the development of circum- 
ferential flow is accounted for in the analysis for stiffness and 
damping. 

One of the problems involved in using and/or understanding 
Black and Jenssen 's results is that various ad hoc governing equa- 
tions are developed and used to obtain separate results. Further/ 
the governing equations do not generally reduce to recognizable 
turbulent lubrication equations. Allaire and Lin [10] have 
partially remedied this deficiency by performing a numerical 
analysis for short seals based on Hirs ' [11] / [12] turbulent lubrica- 
tion theory. However, they have retained Black's initial assump- 
tion of a one half speed circumferential velocity. 



The present analysis also begins with Hirs ' governing equa- 
tion, but has the objective of developing analytical expressions 
for the seal dynamic coefficients incorporating all of Black and 
Jenssen's various developments. The development is based on a 
perturbation analysis of Hirs' equations in the eccentricity 
ratio e, and yields results for a centered zero-eccentricity 
position. 



Governing Equations 

Figure 2 illustrates the seal geometry. Hirs' governing 
equations are provided in Appendix A and are thoroughly discussed 
in references [11} and [12] . We propose to expa^iu these equations 
in terms of the perturbation variables 
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to obtain a "short bearing" solution. The short-bearing aspect 
of the solution is obtained by setting |;^ = 0 in Eq. (A. 2). Sub- 
stitution of the above variables into the governing equations 
(A-l) through (A. 3) yields the following zeroth— order axial 
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The zeroth-order continuity equation is trivially satisfied. 
The corresponding first-order axial momentum equation is 
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while the first-order circumf erential-motnentum equation is 
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and the first-order continuity equation is 
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We begin the solution of the equations by introducing the 
variables 


Uoo = ^ / 2 = z/L / u - = = b 

90 2 zO Roj R^q 


( 7 ) 



into the zeroth-order circumferential momentum Eq. (3) to obtain 
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n - m^ 

®l“2b^C^^C0 


l+mp 

+ b2) 2 ] 

( 8 ) 

(9) 


The variable v defines the "swirl" of the circumferential flow in 
the sense that v(Z) = 0 implies that shear flow is uniformly 
established throughout the seal, i.e., Uqq(Z) = UUqq = Eq. 

(8) defines the development of circumferential velocity as a fluid 
element proceeds axially along the seal. It has the steady-state 
solution v = 0 . 

. For small v, the linearized version of Eq. (8) is 
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By comparison, Black 18] obtains the following definition 
1 tm * 1 


a[l + 8*(— ^)] , 8* = 


l+{~h) ^ 


, mQ = -0 . 25 


( 12 ) 


The parameters 8 and 8* are different because Black retains 
Yamada's 1/7 power law velocity distribution, while Hirs does not. 



However, the two parameters coincide at b = 0 and and are not 

markedly different over the complete range of b. The factor of 

1+m 

two difference between g (1+mQ) and 3*( — follows from differ- 
ences between Hirs * and Black's governing equations. 

As things turn out, the solution to the linearized Eq. (10) 

V = (13) 

is a very good approximation to the decidedly nonlinear Eq. (9) , 

particularly for small magnitudes of v^, the initial swirl. 

Figure 3 illustrates the maximum percentage error in v for Vq “ -0.5 

and L/C = 100 (L/D = 1/2, C/R = 0.01) as a function of b for 

R „ = 3 000 and R = 200,000. The worst error occurs at Z = 1 
au ao 

and was determined by numerically integrating Eg. (9) , and com- 
paring the result to the approximate solution of Eq. (13) . The 
percentage error drops sharply as v^ is moved from -0.5 towards 
zero. Figures 4 illustrates both a and e as a function of b 

for R^« = 3000 and R „ = 200,000, and demonstrate that initial 
aO aO 

swirl will generally persist throughout the seal length. These 
results confirm Black's finding [8] that the 1/2 angular velocity 
assumption for circumferential flow throughout the seal is question- 
able in high axial flow seals if the inlet circumferential velocity 
differs substantially from the assumed 1/2 speed condition. In 
practical terms, neck -ring seals on pump impellers would normally 
have an imposed inlet circumferential velocity close to Ru)/2, 
while interstage seals that are downstream of diffusers would have 
near-zero inlet circumferential velocities. 



Returning to the zeroth-order axial momentum Eq. (2), and 


substituting from Eq. (7) for yields 


= -pa V=^{1 + v^S(mn+l) [1 + - 1)] + 


Eq. (14) defines the nominal pressure gradient between two con- 
centric rotating cylinders and should be comparable to Yamada's 
experimental result. In fact, Hirs [11] makes a comparison to 
Yamada's experimental results assuming = 1/2 (v. = 0). Hirs' 
parameters are n^ = .066. = -0,25 as compared to Yamada's 

values of n^ = .065, m^ = -0.24. Analysis of Yamada's experi- 
mental apparatus and results indicate that the terms in Eg. 

(13) would yield a decrease in calculated values of Xq of approxi- 
mately 2% at R^q - 10,000; R^q = 20,000. For this Reynolds number 
set, Hirs' definition for X^* is approximately 20% too high. Hence, 
including v makes a small, but largely insignificant improvement 
in the correlation. Given that X^ is strongly dependent on R^q 

and relatively independent of these results are not surprising. 

First-order Perturbation Solutions 

In the preceding analysis of the zeroth-order perturbation 

equations, terms of order v^ and higher v;ere shown to be negligible. 

If terms of this order are also neglected in the first-order axial- 

momentum Eq. (4) , the result may be stated 
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Similarly, the circumferential-momentum Eq. (5) reduces to 
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where 


A = Bd+m^) [1 - B(l-mQ)]/[l + Bd+m^)] 


Finally, the continuity Eq. ( 6 ) is stated 
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T'he axial-momentum Eq. (15) is similar to those developed by 
Black and Jenssen [1-4,8]; however, the term which appears 

in the second product on the right-hand side is 2 U^^/b in the above 
references. This discrepancy follows directly from differences 
between the governing equations of Hirs and Black. 

Preceding analyses for seal stiffness and damping coefficients 
[1-4,8,10] have generally* neglected the first-order perturbation 
in the circumferential momentum, which eliminates Ug^ as a variable. 


*Black and Jensen [3] complete an approximate numerical solution 
including a small perturbation about a U 0 „ ^ assumed circum- 
ferential velocity. 2 



and permits integration of Eq. (17) for which may then be 

substituted into Eq. (15) to define -gY"* Eq. (15) is then 
integrated to define which is in turn integrated to determine 
seal reaction forces and stiffness ana damping coefficients. In 
the present circumstances, one would prefer to solve the linear 
variable-coefficient Eqs. (16) and (17) for both Uq^ and 
substitute them into Eq. (15) , and proceed to solve for p^, etc. 
Unfortunately, efforts to analytically accomplish this preferred 
course of action have as yet been unsuccessful. A n\imerical 
solution comparable to that employed by Blaclc and Jenssen [3] or 
a finite-difference procedure is apparently required to account 
for Uo- . In the analysis which follows, U., is dropped and 
Black's initial analysis procedure followed. The general good 
agreement between Black's experimental and theoretical results 
certainly support this approach. 


With Uoi eliminated, Eq. (17) is integrated to obtain 
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Substituting these results into Eq. ( 15 ) yields 
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Integration of this equation yields 
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The boundary conditions to be satisfied by the pressure are 


p. " p(0,e,T) = I u2(0,6,t) (1+C), P(L,0,t) = 0 

1 2 Z 


(23) 


These equations yield the following perturbation- variable boundary 
conditions at the seal inlet and exit, respectively 
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Substituting from the first of Eq. (24) into Eq. (21) , and evaluat- 
ing the resultant pressure definition at Z = 1 yields the follow- 
ing definition for 


at Z 
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Time is not explicitly present on the right-hand side of this 
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equation; hence, the partial derivative U^io 
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arbitrarily be set to zero, which leaves a first order linear 
equation of the form 
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The requirement of continuity, U^^q(9) - U^io^® 2 it) , implies 
-0, and yields the final solution 
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where 


a = ba/(l+^ + Ba) 


Substitution for and Pj^(0,6/t) into Eq. (21) completes the 

definition for p^(Z,9,t). 

Seal Dynamic Coefficients 

The components of the reaction forces acting on the seal 

journal are defined by the integrals 
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The clearance h is defined in terms of the components of the 
seal journal displacement vector by 

h - C~Xcose - Ysine 

Hence 
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etc., for the remaining derivatives. Substituting for (-^) and 
its derivatives into Eq. (27) yields the following definitions 
for the force components 
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If Vq a 0, these equations may be stated in the following format 
comparable to Black and Jenssen [3J 
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Figures 5 illustrates both these functions and the comparable 
functions p*/ p*/ p* from Black and Jenssen [3] . The functions 
Pq/ Pj_ are quite comparable to p*, p£, but are less sensitive 
to changes in b. For b = 0, Pq and p* coincide and the differences 
between p^ and p* are minimum (with respect to b). Although the functions 
yg an(d y| differ significantly in magnitude, particularly for small values o-^ 
a, these parameters generally have a minimal influence on rotordynamic 
response. The fact that P 2 is larger than y^ implies simply that a Hirs- 
based model predicts larger "added-mass" values than those predicted by 
Black and Jenssen. 



Numerical Results 


As noted earlier, experience [5] has demonstrated that inter- 
stage seal characteristics are crucial to the rotordynamic stability 
of the SSME (Space Shuttle Main Engine) HPFTP (High Pressure Fuel 
Turbopump) . The present interstage seal configuration has three 
annular segments at different diameters separated by two seeps. 

A proposed alternative configuration eliminates the steps to yield 
a single-segment seal with the following nominal dimensions 

D = 7.98 X 10~S L - 0.432 x 10“S C = 1.397 X lO"** 

The nominal operating points of the turbopxmp are FPL (Full Power 
Level) , RPL (Rated Power Level) and MPL (Minimum Power Level) with 
the following speeds, density, viscosity, and AP conditions: 
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AP 

FPL 

37,360 

70.28 

1.160 X 10“^ 

1.492 X 10^ 

RPL 

• 

35,020 

68.0 

1.086 X 10”® 

1.314 X 10"'' 

MPL 

23,720 

58.0 

0.759 X 10”® 

7.260 X 10® 


Table 1 contains the seal dynamic coefficients corresponding to 
these conditions for v^ ■ 0 and v^ = -0.5. Inlet swirl causes a 
sharp drop in the cross-coupled stiffness and damping coefficients 
k and c, and a slight increase in the direct stiffness coefficient 
K, with C and m unchanged. Stability predictions for the HPFTP 
would obviously be substantially improved if the coefficients of 
Table 1(a) are employed rather than those of Table 1(b), because 
of the markedly smaller values of cross-coupling coefficient k. 
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Table 1 (a) . 
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HPPTP dynamic coefficients for 
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Table 1(b) 
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Appendix A ; Hirs Turbulent Lubrication Equations 

Hirs* turbulent lubriccition equations [11] , [12] is a bulk-flow 
theory which does not explicitly make any assumptions concerning 
either (a) local flow velocity due to turbulence, or (b) the 
shape of average flow-velocity profiles. Only the bulk-flow 
relative to a surface or v'^all and the corresponding shear stress 
at that surface or wall are considered or correlated. Hirs' axial 
and circumferential momentum equations can be stated, respectively, 
as 
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with the bulk-flow momentum equation 
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